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ABSTRACT 

A technique for calculating ballistic density from Satellite Infrared 
Spectrometer (SIRS) radiance data has been developed and evaluated. The 
experiments with both simulated radiances and real-time radiances show the 
technique has considerable skill in computing ballistic density. 
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ABSTRACT 



A technique for calculating ballistic density from Satellite Infrared 
Spectrometer (SIRS) radiance data is described. A model experiment in which 
radiance values are calculated from the temperature sounding is used to 
verify the technique. In a second experiment real-time SIRS radiances were 
used with synoptic and geographically-consistent soundings for evaluation 
of the ballistic density. For the limited data sample tested, the 
technique of computing ballistic densities from satellite radiances is 
capable of matching those computed from conventional radiosondes. 
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Introduction 



The success of the temperature profile detemination technique (Smith, 
et, al . , 1970) from the Satellite InfaRed Spectrometer (SIRS) data has 
important implications for ballistic calculations. Availability of these 
complete soundings allow calculation of ballistic parameters in remote regions 
without radiosonde soundings. However Commander K. Ruggles noticed an 
analog between the "weighting functions" for the SIRS channels and the 
ballistic density weighting function. The SIRS "weighting function" is the 
logarithmic pressure derivative of the transmissivity (t) of the atmosphere 
for that wavelength. These curves of di/d In p are exact for the partial 
pressure of CO 2 and are related to the total pressure through the constant 
mixing ration of CO 2 . The similarity of the ballistic and transmissivity 
weighting functions suggests that the ballistic density might be determined 
directly from the radiance values. Ruggles (see Appendix A) has developed 
mathematical relations that illustrate the analogy, and give a basis for a 
direct statistical determination of the ballistic density from SIRS radiances. 
The objective of this research is to verify the suggested technique. 

The first check of the technique used 100 soundings (Wark, et. al., 1962: 
Appendix A) . The radiance in each channel which the satellite would detect 
was computed for each temperature sounding using 



I(V ) = B[V T(p )] T (V , p ) 
1 IS IS 



t(v ,p ) 



r 



B[V^, T(p)] dT(v^,p) 



( 1 ) 



with the Planck function 



B[v, T(p)] 




( 2 ) 



-5 2 -1 -1 

where = 1.19064 x 10 erg cm sec (ster) 

= 1.4387 cm °K 

In these expressions v is the reciprocal of the wavelength and T is the 
temperature. Details of the radiance calculations are described in Appendix B. 
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Values of the ballistic density were also computed for each of the 
temperature profiles. The ballistic density was defined 

r^T F (z) dz (3) 

Jo 

T 

F (z) dz 
P 

where Pg(z) is the density distribution of the standard atmosphere at 45N 
for January, ballistic weighting function which was taken 

from Fig. 16b of Finke (1969), and was 66.5 km. Further discussion of 

the ballistic density calculation is given in Appendix C. 

Corresponding values of p and the simulated radiances at the eight 

D 

SIRS wavelengths were statistically screened with the BIMED 02R package 
(Dixon, 1966) using 50 soundings for the dependent sample. The coefficients 
a(v^) in the resulting linear expression 

8 

p = a + 2 a(v ) I(v ) (4) 

i> o . - 1 1 

1=1 

were statistically checked against the data sample p_, I- , ..., Iq 
resulting from the 50 independent soundings. From this first, or simulated, 
approach it was determined that the coefficients in the regression Eq. (4) 
could be obtained with a high degree of stability. 

In view of the successful verification with simulated radiances of the 
mathematical-statistical approach suggested by Ruggles, a second check of 
the technique was undertaken using real-time SIRS radiances in conjunction 
with synoptic and geographically-consistent soundings for evaluation of the 
ballistic density. The general approach was to evaluate the ballistic 
densities with Eq. (3), interpolate to the sub-satellite point, and statistically 
screen as before to obtain a relation as in Eq. (4). The regression co- 
efficients again showed considerable stability. However the so-called 
factor-method of obtaining orthogonalized combinations of the I^, ..., Ig 
as predictors for day-to-day use will also be checked for improved stability 
of the coefficients. 

2. Results with Model Atmospheres 

The test with model atmosphere soundings from Wark, et. al. (1962) is 
indicative of results to be expected with no cloud contamination effects. 
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Temporal and spatial differences between satellite and radiosonde soundings 
are also eliminated by computing the radiances directly from the temperature 

soundings. The transmissivities T for eight channels centered at 899, 750, 

714, 706, 699, 692, 677, 669 cm for use in Eq. (1) and Eq. (2) were provided 

by courtesy of Dr. W. L. Smith of NOAA. 

3 

Application of the stepwise regression procedure to y = (p -1) x 10 

1 B 

as predictand and to the computed radiances for the 50 soundings in the 
dependent set as predictors led to the equation 



y, = 226.26 + .04878X, - 1.29215X^ + .580582X, 

1 13 4 (5) 

- 1.22287X^ + .46124X^ + .07099X^ 

D 7 o 

th 

Here X. is the factor rr times the radiance in the i channel, beginning 
with 899 cm . The multiple correlation coefficient of Eq. (5) is .9978 
indicating the right side of (5) accounts for 99.56 percent of the variance 
of y^* The two variables and X^ were statistically redundant and 

were screened out of (5) . Supplemental stepwise regression equations relating 
X^ to the remainder of the predictors resulted in a multiple correlation 
coefficient of .999. An identical result was obtained using X^ with the 
other predictors in Eq. (5). Thus inclusion of X^ or X^ could provide 
no additional information above that provided by the six-predictor Eq. (5). 

A summary of the statistics for the two sets of model atmospheres is 
given in Table 1. The multiple correlation using the soundings of the 

TABLE 1. Effectiveness of Eq. (5) when used in the dependent and independent 
samples . 





Samples 


Std Dev 


Multiple 


^Expl . 


Std Error 


Max Error 




Mean p 

B 


of Sample 


R 


Variance 


of Estimate 


of Estimate 


Dep 


.98657 


.02380 


.9978 


99.56 


.00168 


.00358 


(n = 50) 














Indep 


.98296 


.02319 


.9942 


98.83 


.00253 


.00689 



(n = 50) 

independent sample to compute the X^ on the right side of Eq. (5) was .9942. 
This means the prediction equation proved to be extremely stable when 
tested against the independent samples, as 98.83 percent of the variance 
was explained. Both rhe standard error and the maximum error of the estimates 
increased in the independent sample. However the values are rather small 
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compared to the standard deviation of the sample. The capability of 
explaining a large fraction of the variation and small standard error in 
the simulation experiment certainly justifies the extension to operational 
SIRS radiances. 

3. Results with Actual SIRS Data 

In this phase of the research, actual values of SIRS-A radiances were 
used rather than computed values. The observed radiance values include 
contributions from clouds in the f ield-of-view of the SIRS spectrometer, 
"Clear-column" radiances obtained by adjustment for cloud contamination were 
provided by Dr. H. Woolf of NESC. Radiosonde soundings (1200 GMT) over 
Eurasia for a two -week period were provided by P. Lowe of NWRF. Eight days 
with corresponding radiance and radiosonde data were available. Unfortunately 
there was a two-day data gap between the four-day periods, December 27-30, 

1969 and January 2-5, 1970. The first period was in the early stages of a 
stratospheric warming that was well-established by the end of the second period. 

Fields of ballistic density computed from Eq. (3) were analyzed from 
1200 GMT soundings, and values were extracted at the locations of the 
satellite radiance data. Radiosonde data lagged some of the SIRS observations 
by up to four hours. However these lags were accepted to provide sufficient 
data for statistical tests. On any one day the number of scan spots, and 
of corresponding p , varied from 15 to 24. Thus for a statistically 

D 

suitable dependent sample, it was necessary to pool the data of the first 
three days of each period. Data from the fourth day was used as the indepen- 
dent test of the regression equation. Because of the data gap, only two 
such independent tests were possible, rather than the optimum of seven tests, 
if no gap existed. 

Another source of unreliability might occur from a bias introduced by 
the correction of the observed radiances for cloud contamination. While 
the details are not important here, it should be noted that the Smith (1970) 
scheme for removing cloud effects uses information from three of the eight 
channels. An iterative scheme is used which determines the elevation of 
cloud layers and their amounts. The procedure minimized the difference 
between the actual radiance and an equivalent dear-column radiance. If 
the iterative scheme suggests nearly overcast conditions at upper levels, 
the sounding is tagged as being less reliable. Nevertheless for the current 
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research all reported equivalent dear-column radiances were accepted as 
valid. 

For both data periods the first three days provided a stepwise-screened 
regression equation, much like Eq. (5). Only five channels were significant 
in each case when the critical F-test of Miller (1962) was used. The pre- 
dictors X^, X^, X^ and Xg provided the best set of five-predictor 

equations, and the best verification when applied to the independent samples 
of data. The results of the tests are summarized in Table 2. 

TABLE 2. Statistics of the three-day dependent prediction equation and 
resulting verification on the fourth day in each period 
(predictors X^, X2, Xg, X^, Xg) . 





Mean 


Std Dev 


Multiple 


Expl 


Std Error 


Max Error 


Period 






R 


Variance 


of Estimate 


of Estimate 


I , Dep 


1.01463 


.00955 


.8832 


78.00 


.00482 


.01022 


(n = 54) 


I . Indep 


1.01821 


.01159 


.9382 


88.02 


.00398 


.00816 


(n = 24) 


II. Dep 


1.01638 


.01616 


.9819 


96.42 


.00318 


.00996 


(n = 69) 


II. Indep 


1.01062 


.01865 


.9636 


92.85 


.00516 


.01107 



(n = 16) 



The best results were obtained in the second period with a multiple 
correlation approaching those found in the simulation experiment (Table 1) . 

The standard error is also increased to .003 - .005, or a percentage error 
in determining the integrated ballistic density of only .3 to .5 percent. 

It should be noted in Table 2 that the dependent regression equation in the 
first period explained only 78 percent of the variance in p„ . However the 

D 

explained variance using this same equation as an estimator in the independent 
test improved to 88 percent. This improvement could have resulted from 
relatively less noise arising from error sources described above. Another 
possible interpretation is that density changes accompanying the stratospheric 
warming were subject to a relatively large temporal variation during the 
three days comprising the dependent sample. 
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4. Conclusions 

The experience with the actual data suggests difficulty in interpolating 
the ballistic density values closer than + .003 with good radiosonde 
soundings. About one half of the European soundings either terminated in 
the troposphere, or were missing data at various levels. This difficulty 
in specifying the dependent data, especially with the time lags involved, 
also affects the results of this research. Nevertheless the indication is 
that the technique of computing ballistic densities from satellite radiances 
is capable of matching those computed from conventional radiosondes. 

Standard errors of the estimates are of the same order as the accuracy of 
analyzing the ballistic density field. 

Some of the errors in the estimate of ballistic density may be attributed 
to the uncertainty in the cloud contamination correction. If such a bias 
does exist, and could be reduced or eliminated, the results with the simplified 
soundings would be indicative of the potential accuracy of the statistical 
technique. The results suggest that as few as five channels might be required 
if the statistical technique was adopted. However with the present instrument 
the remaining three channels are necessary to remove cloud contamination 
effects. 

5. Suggestions for Further Research 

Ruggles also suggested (see development in Appendix A) that a ballistic 
wind vector could be statistically determined using radiances in 

gradient form 



In this equation g is the acceleration of gravity, f is the Coriolis para- 



b^ would be calculated as above for the ballistic density. 

With the present NIMBUS III data spaced nearly 30 degrees in longitude, 
either the ballistic zonal wind component or that normal to the path of the 
satellite could be determined. The requirement for specifying gradients in 
Eq. (6) greatly reduced the number of samples available for a statistical 
test. A more complete test for the ballistic wind is recommended with data 
over North America to ensure sufficient rawinsonde coverage of a high quality. 
It is recommended that "simultaneous" rawinsonde and SIRS-B data be obtained 




B 



( 6 ) 



meter, V is the horizontal del operator, and the new regression coefficients 
H 
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during the transition period from summer stratospheric easterlies to autumn 
westerlies. Availability of rocketsonde data for use in continuing the 
wind profile to upper levels is considered mandatory. It may be necessary 
and advisable to interpolate in both space and time to obtain adequate data 
for use in Eq. (6). A SIRS-B data sample having lateral scan information 
would allow computation of the meridional component as well. 

The expression within parentheses in (6) can be interpreted as an 
"equivalent ballistic geopotential", that is, the integral of the geo- 
potential weighted by the ballistic wind weighting function. Since gradients 
of the geopotential are desired, the actual values of the geopotential may 
be replaced by deviations from the standard heights to reduce the truncation 
error of the integral. Some verification of this "equivalent ballistic D 
value" approach is given in Table 3, which refers again to the simulated 
data. It should be noted that only the D values in the free atmosphere are 
included; the surface pressure is assumed to be the standard value. A 
surface pressure analysis would be necessary to "anchor" the statistical 
estimation of the "ballistic D value". 

TABLE 3. Effectiveness of statistically determining equivalent ballistic 





D values 


(meters) from 


simulated 


radiances . 






Sample 


Mean 

D 


Std Dev 
of D 


Multiple 

R 


io Expl. 
Variance 


Std. Error 
of Estimate 


Max Error 
of Estimate 


Dep 

(n = 50) 


146.47 


336.60 


.9981 


99.62 


22.28 


95.52 


Indep 


209.73 


292.70 


.9978 


99.56 


19.58 


52.98 



(n = 50) 

At least for these simulated radiances the stability of the regression 
coefficients is quite high as more than 99 percent of the variance is 
explained in the independent sample. In addition the standard error of the 
estimate is 20 meters, which is quite low. Of course it is the gradient of 
the ballistic D value which is required in Eq. (6), and the error might be 
larger when adjacent soundings are used. Tests with actual radiance data 
are necessary to check the consistency and the accuracy of the D values. 

In the proposed extension of the research with new data, both the direct 
method with Eq. (6) and the indirect method using the ballistic D value are 
to be evaluated. The latter method is analogous to that presently used to 
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compute ballistic winds from analyzed height fields. It allows a separate 
determination of the capability of determining ballistic wind components 
from a two-dimensional field of the integral of weighted geopotential. 

In another aspect of the research, new regression techniques have been 
evolved, but have not yet been used. These new methods involve factor 
analysis and produce linear combinations of all eight of the dear-column 
radiances as predictors. The predictors thus formed are linearly independent 
of each other. This procedure has not been applied in the current study 
involving , because we felt the standard stepwise regression analysis 

D 

was sufficient vertif ication of the basic hypothesis. However the advantage 
of the factor analysis is that it gives a unique decomposition, and thus 
the coefficients should have greater temporal stability. The availability 
of the new statistical technique and of the new data should allow a more 
complete test of the applicability of satellite radiance data in determining 
ballistic density and wind. 
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Appendix A - Use of Atmospheric Infrared Spectroscopy for 
Ballistic Wind and Density Determination 

after 



CDR K. W. Ruggles (USN) 



The application of atmospheric infrared spectroscopy to the ballistic 
wind and density determination problem is examined from a theoretical basis 
The close relationship between the response function of the infrared satel- 
lite sounder and the atmospheric ballistics is shown. 

Ballistic density for a vertical column of atmosphere can be expressed 



B 

P 



G(Ps) 

In P 



oo 



p(ln p) Wp(ln p) 



d In p 



1-A 



where: 



B 

P 



is the ballistic density of the atmosphere 



G(p ) is a constant based on the ballistic characteristics of a 
s 

standard atmosphere 

W^(ln p) is the missile density weighting factor 
p is pressure 

is the pressure at the surface of the earth. 



In a similar fashion, a ballistic wind for a vertical column of atmos- 
phere can be expressed: 



B 

V 



1 ^ 

In p 



oo 



oo 



jV (In p) 



(In p) 



d In p 



where: 



B^ is the ballistic wind of the atmosphere 
|V is the horizontal vector velocity 



2-A 



9 



As a first approximation, we assume geostrophic flow. Equation 2-A then 
can be expressed: 



oo 



In p 



£ 



oo 



I /k X V^Z(x, y. In p) W^(ln p) d In p, 



3-A 



where g is the acceleration due to gravity, V Z is the horizontal gradient 

n 

of height measured along a constant pressure surface, and f is the coriolis 
parameter. Regrouping and reorganizing, equation 3-A can be expressed: 



= 






oo 



X7 - ^ 1 - Ik X V 

V f In p H 



oo 



Z(x, y. In p) W^(ln p) d In p 



4-A 



Examining equation we shall express the weighting function as 

a linear combination of arbitrary functions, F^(ln p) such that: 

N 

W (In p) = Z a.F.(ln p) . 5-A 

^ i=l 

a^ are constants. At this point the constraints on F^(ln p) are that they 
are functions of In p, and they are known, fixed functions. In a similar 
fashion, we can express p(ln p) as a linear combination of arbitrary func- 
tions, H^(ln p) , such that: 

N 



p(ln p) = Z r H (x, y. In p) . 



i=l 



1 1 



6-A 



where r^ are weighting constants. Substituting into equation 1-A and re- 
arranging, we can write: 



B = 

P In p 



N 
Z 

oo i =1 




oo 



H^(ln p) F^(ln p) d In p. 
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Following similar reasoning, equation 4-A can also be written: 



N 



f In p 



I k X V Z s .b. 

rl • M XX 
OO X=1 



^oo 



y. In p) F^(ln p) d In p. 8-A 



where b^ and s^ are scaling constants related to F^ and respectively. 
Equations 7- A and 8-A are now in a form suitable for examination in compari- 
son with the radiative transfer equation of interest. 
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A simplified form of the radiative transfer equation will be examined 
for one of the spectral intervals of a satellite spectrometer: 

A 

O 



1 . 

1 



^A, f, T(ln p), T(ln p) J + 



B[f, T(ln 



d In p 9-A 



oo 



where: 

is the intensity of radiation sensed by the i-th channel of the 
satellite spectrometer. 

is the contaminant intensity, such as the surface or clouds; 

^ and is a function of A, the amount of contaminant, f, the 

frequency of the i-th channel, and T(ln p) the temperature of 
the contaminant. 

B[f,T(ln p)] is the Planck radiance at frequency f and temperature T. 

T(f, In p) is the fractional transmittance of the atmosphere in the 

spectral interval centered at frequency f, and from pressure 
level p to the satellite. 



Provided the cloud contaminant term can be evaluated and eliminated, equation 
9-A reduces to: 

|0 




B^[f, T(ln p)] 



d In p) 

d In p 



d In p 



10-A 



oo 



Using equation 5-A, we can express a ballistic weighting function as a 
linear combination of spectrometer weighting functions, such that: 



F . — , - 
1 d In p 

Similarly, since B^ is a function of T, and therefore a function of Z or p, 
we can set: 

H. = B. 

1 1 



Equation 10-A then becomes 



I. = 
1 



H.F. d In p. 
11 



oo 



11-A 



11 



Substituting equation 11 -A into equations 7- A and 8- A, we finally show 
that the ballistic density and wind may be directly expressed in terms of 
measured SIRS intensities; through the regression coefficients a^ and b^ , 
respectively, which follow: 

N 

B = G(p ) Z a.r.I- 12-A 

P s 1 i 1 



N 

B - § k X V Z b.s.I, 13-A 

V £ 1.1 1 1 

where N now is the number of radiometer channels used to determine the 
ballistic parameters. 



12 



Appendix B 

Radiance Calculations with Known Model Soundings 



The transmit tances p) for each of the eight SIRS channels are 

shown in the illustration in Table B-1, Transmittances for the seven 

channels in the CO^ band were computed by Dr, W. L, Smith of NESC for 

pressure levels up to .01 mb. For the water vapor window channel, the 

transmittances was assumed equal to 1.0 above 200 mb and calculated from 

200 mb to p = p following Saiedy and Hilleary (1967) 
s 



are shown in Table B-1. 

Temperature values were generally at mandatory reporting levels in the 
library of one hundred model atmospheres drawn from Wark et al (1962) for 
use in the simulation experiment. Any available significant level data was 
also used in linearly interpolating the temperature values to the pressure 
levels with t values (see Table B-1), and the Planck function was calculated 
from Eq. ( 2 ) at each of these levels. The surface pressure (p^) was taken 
to be 1000 mb, although the actual surface pressure was normally higher. 
Simpson's method of numerical integration was used in evaluating the integral 
term of Eq. (1) for the 34 layers between the pressure levels listed in 
Table B-1. 

Although the model atmosphere soundings were carefully selected, few 
soundings reached above 10 mb. Above the final reporting level the tempera- 
tures were extrapolated following the January (45N) standard atmosphere 
lapse rate. This resulted in a uniformally warmer or colder atmosphere, but 
is probably not a serious error in the simulation experiment as both the 
radiances and the ballistic densities are calculated from the same sounding. 




2 -1 

In this equation K = .095 cm gm and q is the mixing ratio of water 
vapor. Transmittance values for this channel and one of the CO^ channels 
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Table B-1 





Transmitt anc e s 


for V = 699 


cm ^ and 899 cm ^ 




Pressure 

(mb) 


''^699 


'’’899 


Pressure 

(mb) 


”^699 


= 1000 


.0004 


.9260 


20 


.9032 


800 


.0012 


.9413 


16 


.9209 


650 


.0037 


.9707 


13 


.9342 


500 


.0176 


.9895 


10 


.9477 


AOO 


.0529 


.9972 


8 


.9567 


300 


.1409 


.9992 


6.5 


.9635 


250 


.2137 


.9998 


5 


.9706 


200 


.3088 


1.0 


4 


.9754 


160 


.4047 




3 


.9803 


130 


.4904 




2.5 


.9829 


100 


.5883 




2 


.9855 


80 


.6600 




1.6 


.9877 


65 


.7170 




1.3 


.9894 


50 


.7764 




1.0 


.9912 


40 


.8174 




0.7 


.9931 


30 


.8597 




0.4 


.9954 


25 


.8813 




0.1 


.9982 








0.07 


.9986 
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Appendix C 

Ballistic Density Calculations 



The form of Eq* (3) for evaluating the ballistic density 



Jo - ^ 

^s 



I 



( 3 ) 



i 



^ F (z) dz 
P 



is essentially the same as that defined in Appendix A, Finke (1969) has 



intervals in the January standard atmosphere. An empirical function was 
used to fit these values and calculate F in layers bounded by mandatory 



mandatory pressures and are normalized with the factor in the denominator of 
Eq. (3). As differences in layer thickness, computed from the hydrostatic 
equation, are included in this weighting factor, the values do not uniformally 
increase or decrease. 

The mean density p in each of these layers is calculated from the 

temperatures (T ,T ) at the mandatory pressure reporting levels ) 

Ad a B 



for each layer are calculated from the same expression using the January 
(45N) Standard Atmosphere temperature values. Missing data at mandatory 
levels was interpolated from significant level data or adjacent mandatory 
level data. As was noted in Appendix B the data were extrapolated using 
the Standard Atmosphere lapse rate above the termination level of the 
sounding. This procedure produced a colder than normal sounding all the 
way to z^ if the last reported temperature was colder than normal. Some 
experimentation to improve estimate of densities above balloon termination 
levels is suggested, although F^ is quite small above 1.0 mb. The 
ballistic densities estimated by regression from the satellite-measured 






where R, is the gas constant for dry air. Standard density values p 
d s 
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radiances will contain information from this layer and could lead to some 
bias. Soundings with surface pressure greater than 1000 mb were truncated 
at 1000 mb, while those with surface pressure less than 850 mb were 
rejected. Between these limits the temperature data were extrapolated to 
1000 mb using the first two temperatures reported. 
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Table C-1 

Normalized Ballistic Density Weighting Function 
Between Mandatory Pressure Levels Used in this Study 



Pressure 


F 


Pressure 


F 


(mb) 


P 


(mb) 


P 


1000 


.0273 


20 


.0224 


850 


.0470 


10 


.0090 


700 


.1368 


7 


.0065 


500 
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